tt a tt St Ht tt HH HH HH Ha 
tt a tt a a tt i 


T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 
T 


HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HH HHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHHHHHHHHHHHHHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 
HHH HHH 


tt tt tt st 4 Ht Ht Ht Ht Ht Ht HH HH 
AAA a 
| LE) 


LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 
LLLLLLLLLLLLLLL 


*eF ILE*#]D**MTHGATANH 


MM MM 
Lal MM 
MMMM 6MMMM 
MMMM MMMM 
aM 60MM OMM 
mma 66MM OMM 
MM MM 
Lal MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
MM MM 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LL 
LLLLLLLLLL 


LLLLLLLLLL 


AAS tt tt st ts 
A444 St tt Ht 


os 
ta fea 


Ge be By Fay Pens Bene Gand Gene fend fewe foe) Qos ewe 
os be 


de be 


(bana Pens beng tena teed Seon beg feo fs Cems imo Demet fn oe 


— be 
tet ee 
bom bee 


tent bat 


HH HH 
HH H 
HH HH 
HH HH 
HH HH 
HHHHHHHHHH 
HHHHHHHHHH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
SSSSSSSS 
SSSSSSSS 
SS 
SS 
SS 
SS 
SSSSSS 
SSSSSS 
SS 
SS 
SS 
SS 
SSSSSSSS 
SSSSSSSS 


ABRABA 


AAARAAA 
AA AA 
Aa AA 
AA AA 
AA AA 
AR AA 
Ak AA 
AAAAAAAAAA 
AAAAAAAAAA 
AA AA 
AA Ab 
AA AA 
AA AA 


44 
“++ 


AS at se tS St HH 
Aaa st SH ts 


HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HHHHHHHHHH 
HHHHHHHHHH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 
HH HH 


8B 9 
MTHSGATANH ; G-floating Precision Hyperbolic Arctan 16-SEP-1984 01:25:48 VAX/VMS Macro V04-00 
Table of contents 


(2) 47 i petat ted serrens Aa bis 
: ; De ve Par 


) ¢} est ores! 9 
) 154 = G-float ng Precision Meee tboltc Arctangent 


os aay ; G-floating Precision ioorbok ie hexeen i ot tat bY 9}: $3: $8 Paya Macro Vv04-00 Page 1 


EP-1 MTHRTL.SRCJ MTHGATANH. MAR; 1 (1) 
»TITLE MTHSGATANH ; G-floating Precision Hyperbolic Arctangent 
. 1DEN /35003/ : File: MTHEATANH, MAR Edit: PDG2002 ~ 


REAAAAAAAAEAAARAEAAAEEARAAARAAAAAAAAAAEAAEAAAARAAAAAAAARAAAAAAAAAAAAARAEAAEAREREAEE 


** ® 
*® ae 
't COPYRIGHT (c) 1978, 1980, 1982, 1984 BY * 
'® DIGITAL EQUIPHENT eGORPORATION, MAYNARD, MASSACHUSETTS. * 
:# ALL RIGHTS R * 

® 
@ THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED « 
:® ONL ACCORD ANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
:® IN ON OF THE ABOVE COPYRIGHT NOTICE. OF TWAR ANY OTHER * 
'« COPIES THEREOF MAY NOT BE PROVIDED OR OT HERUISE MADE AVAILABLE TO ANY * 
;* OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 

® 
** ® 
*® * 
*® ® 
*® & 
*® * 
*® ® 
*® * 
*® * 
*® * 
** ® 


—ZzOoO oo 


TRANSFERRED. 


THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
OORPORAT Ok NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT 


DIGITAL ASSUMES NO ey, 


BI 
SOFTWARE ON EQUIPMENT WHICH I 


ITY FOR THE USE OR RELIABILITY OF ITS 


IBIL 
H IS NOT SUPPLIED BY DIGITAL. 


PARA RAAAAAAAARAAAALASARALARELAASASASAL SAAR A ARAL ARS A ALARA RRR A AAA A ARRAS SDA DO 


fACiLitys MATH LIBRARY 
"ABSTRACT: 


MTHSGATANH returns the enrtest ing precision hyperbolic arctangent of the 
G-floating precision argument. The call is standard call-by-reference. 


VERSION: 2 
HISTORY: 
AUTHOR: 
Peter D Gilbert, 23-Jul-81: Version 2 
MODIFIED BY: 


SOQOooooooooooooooooooooooooooooooooooooooooooe 
MEW 9 ODNAUE WN 9 COONAN EWN 9 OO NAME WA OOOO En 
Be Se Be Oe Oe Oe Se Oe Oe Oe Ge Se Oe Se Se Se Se Se Se Se Se Se Se 


PPS SF AWIWWWIWAIWAININII GIP PNP POPONONPUNP NN 9 PO Os ws Os es 


SOoOSOOCOCSOSOSOSOSOSOSOSOOSOSSSOSOSOSOOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOSOOSSOS: 
SOoOOoOCOCOCOOoOoOOooOOCOSoOoOoCOOCSoSoOSOSoSOOSOOSOSOoOSSoOS 


SOoOQoooQooooooooooooooooocoooooooo 


eee eee ee Se 


o- 


- ~ MT 
spree SPICER GAESEY 'EREEIMG BES ROTEL Bway HOP || 


-SBTTL HISTORY ; Detailed Current Edit History 


Prec 
aile 
4 
4 
; VERSION 1 
; Edit History for Version 01 of MTHSGATANH 
3 5 01-0 Original July 1981 

i ; Edit History for Version 02 of MTHSHATANH 


; 00 Original 
3 i 001 - Ch hange MOVZBL ‘ yg when accessing ip ae -AL OG. V. POG 2-Dec-1981 
; 2 = Change RSB to RET after error exit. PDG 6-Jan- 


—E 9 
ating Precision Hyperbolic Arctan 16-SEP-1984 01:25:4 AX/VMS Ma 4- 
Tl g ; Declarative Part of Modul ~$Eb=1383 94:33:38 LAN TS  SRESRTNGAT ACH MARS 1 ‘oat 3) 


-SBTTL DECLARATIONS ; Declarative Part of Module 


za 
Po 


MTHSGATANH : Gf 
ethos BECLA 


: INCLUDE FILES: MTHJACKET.MAR 


: EXTERNAL SYMBOLS: 
.DSABL GLOBAL 

BINARY CALLS .CONDITIONALS DEF INIT IONS, EXPANSIONS 
“EXTRN " ARGAAT 


IN 
THSK_ [NV 
SSIGNAL 
-EXTRN MTHSSAB_ALOG_V 
~EXTRN MTHS$AB_G_FHT 


: EQUATED SYMBOLS: 


SOooOSCCOSOOOOOCSOOOOCOOOCSOOOOOOOOSOOOSoOSo 
C8 NAUSEA OOOO NAMES WN 0 OONAU EWN 0 OONOU Ew 
. 


— 
pe OR A EN + tt try tet ett PRN RPS A SN Te A e 


DOOOOOOSOOSOSSOSOOSOOSOSOOSOSOOOSOOOSOOOOOOOOOSOOOOOOoOO 


00 
00 
00 
$0 
0 > MACROS: 
3 3 
00 .MACRO OPDEF X, OP, SH 
00 SOPDEF ADDX  “*XOO@SH+OP,R'X,M'X 
00 OPDEF ADDX3  “XO1@SH+OP.R'X.R'X,W'X 
0 [OPDEF SUBX SKOZSH+OP XM! X 
00 OPDEF SUBX3 “XO3@SH+OP.R'X.R'X,W'X 
0 [OPDEF MULX  “*XO4@SH+OP.R'X.M'X 
00 [OPDEF MULX3 “X05 JSH+OP.R'X.R'X,W'X 
00 OPDEF DIVX  “XQG@SH+OP.R'X.M'X 
00 OPDEF DIVX3 “XO7@SH+OP.R'X.R'X,W'X 
00 OPDEF CVTWX “XOD@SH+OP.RW,W'X 
00 OPDEF POLYX  “X15@SH+OP.R'X,RW,AB 
00 SOPDEF MOVX  *X007D,RQ,WQ : mova 
000 OPDEF MOVAX  “X007E.AQ.WL : MOVAQ 
909 “ENDM 
909 OPDEF 6G, <*X40FD> 
09 -OPDEF ADDX “*X06a8+“X40FD,RG.MG 
0 “OPDEF ADDX3 “X01a8+“X40FD.RG.RG.WG 
“OPDEF SUBX  “X02a8+*X40FD.RG.MG 
0 “OPDEF SUBX3 “X03a8+“X40FD.RG.RG.WG 
[OPDEF MULX “X04a8+*X40FD.RG.MG 
“OPDEF MULX3 “X05@8+“X40FD.RG.RG.WG 
[OPDEF DIVX  “X06a8+*X40FD.RG.MG 
“OPDEF DIVX3 “X07a8+*X40FD.RG.RG,WG 
6 “OPDEF CVTWX “xX ba +*X40FD.RW,WG 
0 PDEF POLYX *X : +*X40FD.RG.RW,AB 
0 PDEF MOVX 2400 D,RQ,Wa : MOV 
EF MOVAX “*XO07E.AQ.WL > MOVAQ 


set to exponent 
set to exponent 


~~ 


Suu | aie 4 ae 
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: PSECT DECLARATIONS: 


186 
0000 108 ° .PSECT _MTHSCODE PIC,SHR,LONG,EXE,NOWRT 

: 3 program section for math routines 
1 
112 


: OWN STORAGE: none 
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BECLARATIONS ; Declarative Part of Modul 6-SEP-1984 THRTL. SRCIMTHGATANH. MAR; 1 (4) 1+( 
11 CONSTANTS: 
U3 
000001FC 119 ACMASK = “M<R2,R3,R4,R5,R6,R7, R8> 
118 ; register entry mask and integer 
1} 3; overflow enable 
1 ? i+ The ry in the constants is used to shift the unbiased exponent 
} gi. ; right 4 places so that the rightmost bit is at bit 0. 
154° 
2B00F EF6 e 423FC 125 tNg-t'3 HI: .QUAD aX2BOOF EF Ge Se FC : (Hi 42 bits of ln2)/16 
FIDADSES 47BC3DA : : § LN2_LO: .QUAD “XF1IDAD5E44/7BC3DA : (Low bits of ln2)/16 
1 1 8 LOGTAB1: 3; Constants foe e gts). Generated using 
a : 3 : eq; , sy 10 © art et. al. (sin(2a) 
3 = 
sgonesy eueoerys Guie 19, ay ceagomssraneonens FCP Sg. tUueeeenmnanee 
. D4 = "Vv. 
2996 138 28 ee 5 1 j - QUAD “HC9P9 G89 ¢23tE ; 6 = 1238869039 335909 
BBAC948 PEER. FE 0 134 -QUAD “XB atgee 55 are 3 C5 = -0.16666645767642529 
Grados Gaoeerey Shab 182 “AUB CHRRAEEROLSSparey | EEE EcEOSSSSMRRINC 
4155555 ORES 4001 $ [QUAD = *X erie ttt aa : C2 = 0. 8 $8 % 1555 
FFOOFFFFE FFFFBFFF 04 138 QUAD “XFFOOFFFFFFFFBFFF : Cl = -0.49999999999999112 
00000000 4 9444 05 139 -QUAD * xo090000000000000 : CO = 0.00000000000000000 
00000 be 129 LOGLEN1 = .-LOGTAB1/8 - 1 3 no. of floating point entries 
05 14 
008 128 LOGTAB: : ; Constants for p(z*z). Generated using 
05 164 3: eq. 6.3.11 of Hart et. al. (sin(ga) = 
05 145 $ 3 - 1)/(b_+ 1) where b = 2%*#(1/7)) 
117401 $F 8 FE7 005 138 -QUAD *X VI ZRO1DGE 1G FE7 3; (5 = 5 ES ree ree ee 4 
A587e0 A73FE 06 14 QUAD “xXOBA58 £8 1A73FEC :¢& = OQ. 18457493 oge¢ 
f 0898 ones re 6 148 -QUAD *xC30898 anse ae: $ C3 = Q. {eee lgce 6551 
43444 | 9 FF 0 149 -QUAD “*X 3444944 O3FF : C2 = 0.399999999996 ree rt 
5560555 Rozne0us 07 150 -QUAD = *X 60555 eT SE : C1 = 0.666666666666667851 
00000000 00004020 008 151 .QUAD Fe eect 0 : CO = 2.000000000000000000 
00000005 008 152 LOGLEN2 = .-LOGTAB2/8 - 1 


Ritstiniiieeaniedih rc 


MTHSGATANH 
thos 


arisen 


00000004 


OOCCGCOCOOOOOOOooO 
‘G2 C0 Cd Co CD CD CD CD CD Cd OD C9 09 CD CD CD 000D 


CDOOOCOOOOOOOOOOOOSOOSOOSOOSOOCOOO 


OOSooOooooooooooa 


oooocsoo 


OOOOoooooooo 


0.000909 00 09 C9 0d C9 00 G9 Co G0.G0 Go On Cd 0d 09 Od Cd Od Cd Od Cd Gd CD CD Gd CD CD CD. 


—DBDOODODOOOS000 000 OOOO 0909 09 69 09 09 09 09 09 09 SIN NI NINN NI NNN AOAAAOAOAAAOOUIINIV «=—-=O 


SONA ALAN 9S ODNA UE WIN SO ODNAU EWN 3 O OO NOAUE WIN 20 OD NOAUE WHO OONOUE 


re 1 Hyper olic aeeten 16-$ 1: AX/VMS Macro V04-00 Page 
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pa MTHSGATANH - G-floating Precision Hyperbolic Arctangent 


FUNCTIONAL DESCRIPTION: 

GATANH = G-floating precision floating point function 

GATANH(X) is computed using the following approximation technique: 
If {Xt >= 1.0, error. Otherwise 
Let (1#X)/(1-x) = f ® (2e#n), where 1/2 <= f < 1 


If n is greater than or equal to 1 then 
et set N=n-=- 1 and Fl = 2ef, 
se 


set N =n and fF = f. 


f iF = 11 < Qeeeh t 
Deatanh(X) = netnc2) + W + WeP(W), 
were We UCLSED/ (1=P))#oeeN - 
and P is a polynomial of degree “tes, D=9. 


se 
Obtain FHI croubty s equal to F) from table Look 
2tatanh(xX) = bnC(14xX)/(1-x)) = Neln(2) + in(FHI) + 2#0(2*2), 
Qisa pot ynonial of Fat F=2,D=5, 
Z = (F = + FHI) 
F = (2@®—-N) #(14X)/(1-xX) 


where 


woke Y = FHI *2*eN 
(14X 
+ Led /UC1+Y 
NOTE: The quantities Ln(A=FHI) and In 
equations in two parts - a high eert (containing the 
h h order bits) and a low part (containing the low 
order bits. In the code the h 9" and Hy parts of the 
Lo gag ere indicated by a _HI and suffix eyages 
tive ly. e values were chosen such thet N*LN2_H 
LN_ FH AHI 7 exactly representable. 
CALLING SEQUENCE: 
gatanh.wg.v = MTHSGATANH(x.rg.r) 
INPUT PARAMETERS: 
X=4 ; Contents of x is the argument 
; IMPLICIT INPUTS: none 


: OUTPUT PARAMETERS: 


in a 
1+X) 
+ X 
+ Y 
+Y) 


NN x<x<~<~<5 


used in the above 


Sete Ge Ge Ge Ge Ge Se Ge Se Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Ge Se Ge Ge Ge Ge Ge Fe Ge Ge Ge Ge Fe Ge Gs Ge Se Ge Ge Ss Ge Ge Ge Se Se Se Ge Ge Ge 
+ 

= 

— 

@ 

i 

@ 
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VALUE: G=-floating precision hyperbolic arctangent of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 
SIDE EFFECTS: 


Signals: MTHSK_INVARGMAT if {Xi >= 1.0 with reserved operand in RO (copied to 
the signal mechanism vector CHFSL_MCH_RO/R1 by LIBSSIGNAL). 

Associated message _ "Floating Overflow in math Library’’. Result is 
reserved operand -0.0 unless a user supplied (or any) error handler changes 


NOTE: This procedure disables floating point underflow and integer 
overflow, causes no floating overflow or other arithmetic traps, and 
preserves enables across the call. 


Note: This routine is written to avoid causing ony integer overflows, 
floating overflows, or floating underflows or divide by 0 conditions, 
whether enabled or not. 


‘G9 G9 OD CD Cd Cd CD CO CD CO CD CD CD CUCD CD CDODODCDODODOCDCDCDCD ) 06 a& 
00 C0 CO. On Ob CP GD CD 0D CD 00 C9 CD. 0D CD CD. CD CD CD CD CD CD CODCOD CD 


Mm Be Se Se Se Se Ge Se Ge Ge Ge Ge Se Se Ge Se Se Sse Ses Sse Ses Sete 


DWONAVUES WN (OOO NAUS UO OONOU ES UT 


WAI WII IAI ROPONIPOPONONININ 3 SS SSS 


0 
O16A 31 RR: BRW ERROR 
O1FC OBB ENTRY MTHSGATANH, ACMASK :; standard call-by-reference entry 
BBD 3; disable DV (and FU), enable IV 
50 04 8C 7D sd re MOVX @X(AP), RO ; RO = arg 
20 7 091 rt: MOVX RO, R6 
52 08 6 43FD 94 4 SUBX3 R6, $*#1.0, R2 3 R2 = 1-K 
Fe 15 99 44 LEQ 3; ATANH(X) is not defined for X>=1 
56 : <or? 98 45 ADOX $*#1.0, R6 
1 F rt BLEQ ERR 3 ATANH(X) is not defined for X<=-1 
2 2 $360 1 4 CVTGF R2, R2 
6 § FD 3 $8 CVTGF Rg. R6 
6 46 4 DIVF2 R2, R6 ; R6 = approximation to (1+xX)/(1-x) 
58 00000000'GF 3 g 0 MOVAB G“MTHSSAB_ALOG_V, R8 
58 s ' ¢ 1 ADDL (R8), RB ; R8 = address of ALOG table 
55 56 OO7F 8F 6 § BICW #iaf frp) R6, RS ; RS = Biased exponent 
55 4000 8F A SuBw 86s #*x4000, R5 : R5 = Unbiased exponent 
6 1 : BLEQ NE P ; Branch to processing for n=<0 
55 080 8F Ag § SUBW #iaF EXP, R5 ; Exponent is positive, RS =N=n- 1 
6 5 C SUBW RS, R6 : RO =F = 
6 6 A cB 8 MOVZBL R6, ; R6 = index into ALOG table 
090 3 CE 59 IF ne F_EXP=-X_EXP 
55 8 cf $9 bivue 7] &<F“EXP-X"EXP>, RS ; Shift RS to scale X-floating 
7E 55 4DFD D1 ¢ CvTWx RS, =-(SP) ; Push N onto the stack 
58 6846 9. DS «626 CVTBL (RB&)CR6), RB : RB = offset into FHI tables 
D9 64 BLSS LN_1 44 ww ; Branch to handle F close to 1 
58 eoccsene "er 3 7E 8 65 MOVAX G*ATASSAB-G_FHICR8], R8 ; RB = Address of FHI 
6 «68 7D + 96 MOVX (R8)+, RO ; R6 = FHI 


MTHSGATANH : G-fl ting Pr ision rbo te” aren 16-SEP-1984 01: AX/VMS v04-00 mT 
itt aris $GA = G- “floating "Thecte sion Hyperb 6-SEP- 1382 9 33; § YRTARTL L.S SRTH HGATANH.MAR; 1 _ &) Ta 
3; Compute Z = (F = of Ht)! (F + FHI) 
; os hs Hy = ¥*(1 mat Shee: + reer, 
3 - Y + x*V¥Jj/(1 
: where Y = Fule2een, roughly equal to (1900/ 
5 5 a0 , ADDW » R6 ; RO = FHI * 2eeN = SFHI 
2 i aa in dr a 8 eo 
Ve e = 
46FD DIVX R6, 3 RS = (1-SFHI)/(C1+SFHI) = D 
56 41FD ADDX3 RO, R2, R6 3; R6 =D + X 
ius BS. a tans 
g 46FD R2, RO t R6 = (D+X)/(14D"X) = 2 


Compute Z**2, P(Z*#2) and Z*P(Z**2) 


50 36 
FF4S CF F3 


aon 
une 
Su 
es Dee es | 
ooVo 


RG. R6, RO : RO = J**2 
POLYX RO, #LOGLEN2, LOGTAB2 ; RO = P(Z2*#2) 
MULX R6,. R + RO = Z*P(Z*#2) 
: Compute B = N*LN2_LO + LN_FHI_LO + Z2#P(2*2) 


MULX3 (SP), LNg.LO, R2 ; Rg N*LN 
(R8)+, R ; 


FW 9 OONA UE WN - OOONOUE Wh "OO 
° 
o 
_ 
< 
= 


OOOO O00 00090909 C99 O98 SII NNO OO 


ee ae ed a a a ee ee ae ee ee ee ee a ed ed ad a od od od od ad od od = I IOV OOOOOOOCOOOCOCO +0 


DPD PPP B BB BS BE AIAN III 0 8 8 9 PQQ 1M 78 $08 SHU FTUETUETUTICTICVICTIOTY +4 


SNF NNN NNN EO BD BPN MMM OOO BAA ANT NINN NII INIT DPM OAOAOAAAO 


SOGCOOOOOSOCOOOCOOOOSOOOCOOSOOSOOSOOSOSOOOSOOSCOOOOOSOOOOOOOOOOOoOOoOO 


* LO 
ADDX : Ro = N*LN -t0 + LN_FHI_LO 
ADDX R2, RO > RO=8B 
34 > Compute A = N*LN2_HI + LN_FHI_HI and ALOG(X) 
52. FED4 cf BE 45FD 36 7 MULX3  (SP)+, LN2_HI, R2 ; R2 = N*LN2_HI 
; § 4OFD 9 ADDX (RB), R : RO = A = N®LN2_HI + LN_FHI_HI 
5 4OFD 98 ADDX R2, RO : RO = A+B = ACOG(X) 
50 10 Ae 99 SUBW2  #1aX_EXP, RO : Divide by 2 
04 et RET 
§ NEG_EXP: 
56 55 A 0 SUBW aR, -R6 > R6 =F = 2¢ 
56 56 9A 04 MOVZBL R6, R6 : R6 = index into ALOG table 
09900005 05 IF NE, F_EXP=X_EXP 
55 0 6 Divwe #IS<FIEXP=X"EXP>, RS = ;_- Shift RS to scale ne 
7E 55 4DED 08 CVTWX RS, =(SP) : Push N onto the st 
58 6846 «(98 99 me mee (RB) CRO], RB : RB = offset into FMI phen 
5 19 11 BCss LN_1 PLUS : Branch to handle F close to 1 
58 00000090"6r 48 7E \¢ MOVAX  G*ATASSAB_G_FHICR8], RB ; RB = Address of FHI 
56 68 7D 13 MOVX (RB), R6 3 Re = FHI 
15 : Compute Z = (F = FHI)/(F + FHI) 
1 : l= C1 sx) - Yr (TeX VET 4X) + ye(1ox)] 
3 Z=(Ci+*+x-Y+ + 
18 : where Y = FHI*2**N, roughly equal to Ce 
6 5 a0 a7 ADDW = RS,_-R 3 RO = FHI * 2e*N = SFHI 
52 6 43FD 1 SUBX3 RO, $*#1.0, R2 : R2 = 1 = SFHI 
40FD : ADDX  $*81.0, R6 + R6 = 1+ SFHI 
2 6 46FD DIVX RG, R : R2 = (1-SFHI)/(C1+SFHI) = 
56 41FD 4 ADDX3 RO, Re, R6 > R6=D +X 
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1 ADbX 8.48, R2 : Re . ; : Dex 
DIVX R2, R6 : R6 = (D4X)/(14D"xX) = 2 


Compute Py P(Z**2) and Z2*P(Z**2) 


50 MULX 3 R Z**2 
FEDS CF 3 se nO. ; MboLEN2, LOGTAB2 n9 = P(Z*#2) 


at 
GATAN 
44FD 01 5 
40FD 01 6 $ 
46FD 
: 5 i 
6 45FD i ? : 
30 2fe 153 g 
6 44FD 4 ; MULX = 2*P(2*#2) 
a Compute B = ia + LN_FHI_LO + Z*P(Z*2) 
52 FETA cf $ 45FD 01 § : MULKX3 (SP) LNZ_LO, R2 ; R2 = N®LN2_LO 
8 40FD 18F 3 ADDX  =(RBS, R : R2 = N®LN2-LO + LN_FHI_LO 
50 «552 «40FD 193 39 ADDX R2, RO > RO=8 
19 41 Compute A = N*LN2_HI + LN_FHI_HI and ALOG(X) 
52 FE63 CF f 45FD 0197 rk . MULX3 (SP)+, LN2_HI, R2 : R2 = N*LN2_HI 
28 7 aged he 44 SUBX =-(R8), R2 : R2 = A ° N®LN2 Ay : tne FHI_HI 
5 32 40FD OIA 45 ADDX R2, R : RO = A+B = ALOG 
50 10 A2 1Ag 46 SUBW2 #1aX_EXP, RO ; Divide by 2 
04 a 44 RET 
TAA 349 ; 
be 2 : Special logic for F close to 1 
1AA 352° 
1AA 2g LN_1_PLUS: 
56 08 50 43FD O1AA 54 SUBX3 RO, $*#1.0, R6 3 R6 = 1-K 
55 8 1AF 55 TSTW RS ; Determine which way to calculate W 
12 #1 181 36 BEQL 10s 
56 10 56 47FD 0183 35 DIVX3 sae. 0, R6 : R6 = 2/(1=x) 
56 8 42FD 0188 358 SUBX Bo5,5 0, > R6 = (14K)/ (1x) 
56 A2 018C 59 SUBW RS : Scale R6 
56 08 42rD 18 60 SUuBK Sos .0, R6 + R6 = 
56 «550 38 47FD 18 ¢ 10$ DIVX3 R6 spre > R6=XK/ (1x) 
56 10 AO O1CA 36 ADDW Sraxte R6 | RO = W = 2eK/(1-K) = (14K) 12K) = 1 
FE3C CF 9 36 SSFD O1CD 4 208: POLYX R6, Pri OGLENT, LOGTAB1 ; RO = Q(W) 
6 44FD 0104 36 MULX R6, R : RO = Wea(w) 
52. FE2A CF : 45FD 0108 66 MULX3 (SP), LN2_LO, R2 > R2 = N®LN2_LO 
4OFD O1DF 36 ADDX R2, RO > RO = N®LNO-LO + WeQ(w) 
28 6 torp 15 g ADDX 6, RO : RO = N®LN as’ + WeQ(W) + W 
3 fe 14 cF 44FD O1E é MULX LN@_HI, (SP) ; (SP) = N*LR2_HI 
3 3 40FD O1ED 0 ADDX + : RO = ALOGEXS” 
0 1 Az IF 4 suBW2 #iax_ExP, RO : Divide by 2 
fs 78 ; 
7 ge 3; X <= 0.0, signal error 
7E 00° BF 9A OTF ds ERROR: MOVZBL #MTHSK_INVARGMAT, -(SP) ; condition value 
50 01 F 79 O1F9 7 ASHQ #15, #T, RO : = resu . = reserved operand -0.0 
1FD oe 3 aes to signal mechanism vector 
1FD 7 3; UCHFSL nc R /R1) so error handler 
1FD $0 3 can modify the result. 
00000000'GF 01 FB O1FD 1 CALLS #1, G*MTHSSSIGNAL 3; signal error and use real user's PC 
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RET 3: re ¥ 
; CHFSL_MCH_RO/R1 


a 9 
MTHSGATANH : G-floating Precision Hyperbolic Arctan 16-SEP-1984 01:25:4 AX/VMS Macro Vv04-00 Page 
Symbol table ° 7 pate 94:33:38 YRTHRTL. SREJATHGATANH.MAR; 1 . 
ACMASK = 
ERR 8 
ERROR 1 
F_EXP = 
LR2_H1 1 
LN2_LO 1 
LN_T_PLUS 1 
LN71 PLUS. 1 
LEN = 

LOGLEN s 

OGTAB 01 
LOGTAB2 5 1 
MTIHSSAB_ALOG _V eeneeee x 
MTHSS$AB_G FHT reeeenee x 
MTHSSSIGNAL eereeree x 0 
MTHSGATANH 00000088 RG 1 
MTHSK_INVARGMAT tteeneee x 0 
NEG_EXP 00000139 R 1 
Xx = 00000004 
X_EXP = 00000004 

wer rose nermoeccese + 
; Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
- «ABS. 00000000 g.} 00 ¢ 0.) NOPIC USR CON ABS '_ LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
- . J VEC L 
MTHSCODE 00000205 ¢ 517.) O1¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
dow meer nee merece n seo ewad 
! Performance indicators : 

ase age faults me apse me 
Ph P faul CPU Ti El d Ti 

nitialization :00:00. 7:00:00. 
Initializati 29 00:00:00.07 00:00:00.85 
Command processing 114 0:00:00.71 38 3893 4.64 
Syabol table sort “ 8: 6°84 6: 8: 4°64 
bees 81 0: 39 (b0s0 108-92 
Symbol table output :00: 0 0:00:00.1 
coset synopsis outeut 0: f : i 0:00:00.02 

ross-reference outpu :00:00. :00:00. 
Assembler run totals 33 8 : 98 80; §:99-88 


the working set Limit was 900 pages. 
7026 bytes (14 pages) of virtual memory were used to buffer the intermediate code. 
There were 10 pages of symbol table space allocated to hold 20 non-local and 2 local symbols. 
445 source Lines were read in Pass 1, producing 11 object records in Pass 2. 
pages of virtual memory were used to define ¢ macros. 
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04-00 Page 1g 
VAK=-11 Macro Run Statistics (5) 
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MTHRTL. SRC JMTHGATANH.MAR; 1 
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! Macro Library statistics ! 


peremecene ces enscsceeomeonawena $ 


Macro Library name Macros defined 
_$255SDUA28: (SYSLIBISTARLET.MLB; 2 Sa nti AF 
0 GETS were required to define 0 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHGATANH/OBJ=OBJ$:MTHGATANH MSRCS$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) +MS 
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